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I, IHTRODDCTION 


The work describt'v ; ierei*> is in support of a contract called "Beneficial 
Uses of Space" which ♦•h** R«»reral Electric Coinpany is carrying out for the 
National Aeronautic:-- j-vce Administration At this tine. 

The objective of the study described in this report is to develop quanti- 
tative relations between various force fields and their effects on the 
notion of particles of various sizes and physical characteristics* More 
particularly, the sti'y is concerned primarily with weak force fields 
whose effects are ordinarily masked out by the large gravitational force 
field in the earthbound enviroraaent. 

In a satellite, howe", where the gravitational field is nullified by 
centrifugal force, these ir> latlvely weak forces are strong enough to pro- 
duce substantial particle .'motions which should be useful in sorting or 
separating or otherwise manipulating a variety of particle types* 

If the particles to be manipulated are floating in an evacuated chamber, 
their acceleration will be the important parameter, -and we have found 
that very substantial accelerations can be achieved in many cases* On 
the other hand, if the particles are confined in a container fil^d with 
a viscous medium such as air or water, they will ordinarily reach a drag 
limited constant velocity quite quickly in response to the applied force 
field, so that velocity rather than initial acceleration is the signifi- 
cant parameter. 

As shown in the body of the report, the particle accelerations produced 
in a vacuum or the particle velocities produced in a viscous medium are 
dependent on the nature and intensity of the force field, the viscous 
medium, if used, and the characteristics of the particles themselves. 
These dependencies offer a variety of potential means of sorting or 
separating mixed assemblages of particles. Combinations may be chosen 
%diich will produce motions primarily sensitive either to electrical 
properties, magnetic properties, size, or combinations of the above 
parameters with density, depending on the nature of the force field and 
medium used* 
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It is hoped that the general expressions developed for these particle 
force field interactions and the numerical examples worked out for typi- 
cal cases will enable potential users of the space environment to iden- 
tify fruitful applications in the manufacture or processing of various 
commodities* 

It was not a part of the study reported here to seek out the applications, 
but rather to produce a "box of tools" useful in the search for applica- 
tions* 

The phenomena which have been explored for use in manipulating particles 
are as follows: 
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• LIGHT 

• HEAT 

• MICROWAVES 

• SOUND PRESSURE 

• PARTICULATE RADIATION 

• ELECTROSTATIC FORCES 

• MAGNETIC FORCES 


Brief Description and Coiapartson of Forces 

It is seen that each of the first three phenomena is a form of electro- 
magnetic radiation . The net DC force produced on particles in the field 
of an electromagnetic wave results from 'the radiation pressure of the 
wave. These three forms of radiation, plus X-rays and gamma rays, which 
are also forms of electromagnetic radiation, are treated together in a 
. single section on electromagnetic radiation effects. The force produced 
is proportional to the effective area of the particles times the inten- 
sity of the radiation. 

Unlike electromagnetic waves, sound waves require a medium for tx^ansmis- 
sion. Like electromagnetic waves, sound waves also produce a radiation 
pressure force which is proportional to the effective area of the parti- 
cle and the intensity of the wave. In most cases, however, the particle 
is much smaller than the wavelength, and in this regime the effective 
size of the particle as far as scattering of the wave is concerned is 
much smaller than its actual size, so that very little force is produced. 
This is discussed further in the section on sound pressure, where another 
effect known as acoustic streaming, which may have more useful applica- 
tions in particle manipulation is also discussed. 

In the section on particulate radiation , the case of an electron beam im- 
• pinging on particles is discussed. As in electromagnetic pressure and 
sound pressure, the force produced is proportional to the particle area 
and to parameters related to the beam Intensity. Electron beams differ, 
■however, from electromagnetic waves and sound waves in that the effective 
cross section area of the particle will generally be approximately equal 
to its actual cross section area as far as varlaticns Ln particle size 
are concerned, so that the equations of motion would be expected to hold 
down to very small particles. 

The report concludes with discussions of magnetic forces and electrostat ic 
forces . We are all familiar with the fact that large forces are developed 
on soft iron parts in relays and solenoids by relatively small electromag- 
nets, but it is not so commonly realized that nearly all substances have 
magnetic permeabilities differing sufficiently from unity to produce appre- 
ciable motions in a magnetic field in a zero gravity environment. The 
equations for magnetic forces and their resulting accelerations in a vacuum 
or velocities in viscous media are given, as well as numerical examples for 
various substances. 
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The electrostatic force of dielectrophoresis is shown to be exactly sini- 
lar in form to the magnetic force. It depends on the strength of the 
electric field and its degree of non-uniformity, the dielectric constant 
of the particle, and the volume of the particle, in exactly the same 
functional way that the magnetic force depends on the strength and the ^ 
non-uniformity of the magnetic field, the magnetic permeability of the 
particle, and the volume of the particle. 

These electrostatic and magnetic forces are dependent on particle volume 
rather than on cross sectional area as in the case of the forces arising 
from electromagnetic waves, sound waves and particulate beams. They, 
therefore, furnish a sensitivity to different particle parameters, a fact 
which adds versatility to the application possibilities. 

Throu;^hout this report we refer to the motion of particles of micron size 
in an evacuated or gas filled chamber resulting from various applied 
forces. One must realize that this type of interaction would in general 
be totally swamped out in the earth's gravitational field, since the 
particles would all fall to the bottom of the chamber, particularly in * 
the case of a vacuum chamber. In the zero gravity environment, however, 
the particles would be floating lazily about the chamber, and the applied 
force would be effective in sorting or other manipulative processes. 

For the case of particles such as viruses which must be preserved in a water 
solution, the effects of the forces studied here are very small indeed because 
of the drag force exerted on the particles by the viscosity of the liquid 
medium. For example, the sumnjary chart on page 39 lists limiting velocities 
for particles in air as the viscous medium for various force fields. Limiting 
velocities ranging from 10“^ to 3 x 10“3 cta/sec are shown for typical forces 
and particle sizes in air at atmospheric pressure and temperature. But in a 
water solution which has viscosity fifty times greater than air, these 
velocities are reduced by fifty fold. In electrophoresis systems, separation 
times of tens of minutes in liquid colloid solutions are not uncommon, 
corresponding to velocities of the order of 10’^ cm/sec, but the predicted 
values for the velocities in liquid solutions for the forces described herein 
would be -generally much lower. 
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II. TirER>tAL MOTION AND VrCSCOUS FORCES 

Before actually d . . Ibfng the details of the five categories of small forces 
that were considev^u la the study, it is worthwhile to discuss the environ- 
mental constraints that will apply in general to all of the forces. The tv:o 
Issues to be conr’i’r::. . ’ 1) the competing effects of thermal motion and 

(2) the viscous f>'>:;i:.'-s that will be acting on the particles if they arc in 
fluid. A particl* - .n’-.-. jte temperature T will have a thermal kinetio 

energy^-t(dv^/^ -ikT, • Boltzmann's constant, M is the particle mass, 

V is the particle velocity, and the brackets indicate an average. If a gas 
of these particles is in a container and a force F is applied in a direction 
perpendicular to a wall of the box, the equilibrium distribution of the 
particles is proporticna’ ,to{xp(if^//tJ|)^here x is the distance away from the 

wall. Approximately /0% of tl»e particles will be located a distance 
away from the wall, much like the fact that 70^1 of the earth's atoosphere 
Is located within 10 km vf the earth's surface due to the force of gravity. 

It should be noted ti:." ; ti.s is an equilibrium distribution (the Maxwell- 
Boltzmann distributici.) and does not contain any information as to the time 
necessary to establish the equilibrium. This time might, in fact, be very 
long if the medium were a viscous fluid. In Figure II-l, the initial accelera- 
tion needed at 300°K to establish a distribution withA'rs /c*»\ is plotted 
(curve (1) ) as a function of particle size for spherical particles with den- 
sities of one gram/cm3. The ordinate of Fig. II-l is the initial acceleration 
that a particle would experience in a viscous fluid as a result of the applied 
force, and is the same as the particle acceleration in a vacuum, F/M* where M 
Is the particle mass. It is used here because much of the subsequent work is 
given in terms of the acceleration of particles in a vacuum rather than the 
force on the particles. As can be seen, the necessary initial ac.cel oration is 
a function of the reciprocal cf the particle radius cubed. 

In many cases the particles of interest will be imbedded in a viscous fluid or 
gas and it is then likely that the need for a reasonable particle velocity 
will determine the necessary force rather than the thermal considerations 
given above. This can also be seen in Fig. II-l where Stokes law*- is used 
to Indicate, as a function of particle radius, the Initial acceleration needed 
fo achieve a final velocity of 0.1 cm/sec for spherical particles of density 
1 gram/cm3 at 30QOK for an air medium at standard pressure (curve (2) ) and 
for a water medium (curve (3) ) . As can be seen, both of these curves indicate 
a much higher needed initial acceleration than given from the thei nal considera- 
tions. Also, as can be seen, this requirement on initial acceleration varies 
only as the reciprocal of the radius squared. These curves are meant to serve 
only as an illustration of the type of dependencies to be expected. For a 
specific application, the curves may require scaling since both will depend oa 
temperature and curve (2) depends strongly on pressure. 

It should be noted that the restriction of particle velocity in a viscous fl-cid 
is not necessarily a detrimental effect. The dependence of drift velocity ia 
a viscous fluid on particle radius and force, F (see Fig. II-l for 
explanation) may actually become an additional means for achieving the sepa- 
ration of particles of different sizes if forces of sufficient magnitude caa 
be achieved. 
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III. RADIATION PRESSURE 

In this section^ we will consider tlie use of radiation pressure as a means 
for moving small particles and also as a means for achieving selective 
separation of small particles on the basis of size of the particles, refrScitive 
index of the particles^ and the wavelength of the electromagnetic radiation. 

We will first derive an equation for the acceleration in a vacuum produced 
by radiation pressure. We will then consider the entire range of electromagnetic 
radiation from microwaves to x-rays and gamma rays and attempt to characterize 
the size of achievable acceleration as a function of wavelength in order to 
determine the most promising sources of electromagnetic radiation for radiation 
pressure. Anticipating the result of this study, we will focus attention on 
the radiation pressure produced by sources in the near ultraviolet, visible, 
and infrared wavelength range. Resonant excitation of atomic absorption 
lines in atoms and molecules will also be treated. 

A. General 

The radiation pressure of an electromagnetic wave on a body that is perfectly 
absorbing is simply W/c,(3/vhere W is the power carried by the wave per unit area 
and c is the speed of light. Expressed another way, W/c is the rate^ that 
momentum can be removed from the wave per unit area of the wave. Thfe force 
on a body is given, in general, as 


F A 


iii-i 


where is the efficiency factor for radiation pressure and A is the cross 

sectional area of the body. Consequently^ Qf>r’A is defined as the cross- 
sectional area of the body that is effective in receiving momentum from the 
electromagnetic wave. C3.r contains the "physics" of the momentum transfer . 
process and, in general, is a complicated function of the parameter 

and the refractive index m of the body, where is related to the size of the body 
and ^ is the wavelength of the incoming radiation. (^r can be calculated 
from the efficiency factors for absorption «nd for scattering 

and is given as 




III-2 


where Q is the angle between the scattered light and the transmitted light. 
This equation indicates that all the radiation that is absorbed contributes 
directly to momentum transfer, but scattered radiation only contributes if it 
is scattered away from the forward direction. For exjnple, perfectly absorbing 
particles that are large compared to the wavelength of the Incoming radiation 
(the geometrical optics requirement) have an effective cross sectional area 
equal to their geometrical cross section This is analogous to an 

inelastic collision in mechanics. A large perfectly reflecting particle has 
Q , 80 that ‘the effective cross section for momentum transfer is twice 
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the geometrical cro:.- section, vhich.is analogous to an clastic collision 
in mechanics. 

The initial acceler*’-'; • o;... • <,celeration in a vacuum for spherical particles 
of ra<Jius a (cm) and density ^ (grams/cm3) in a radiation field of intensity 
W (vatts/cra2) is give ‘ ‘ . 


n F , 


=. id" ctvi in- 3 

o f 


vhere M is the mass of the particle. Although the details of the subsequent 
work are for spherical particles since they are easier to handle analytically, 
the results can be exp*'^i ifu to approximately apply to particles of irregular 
shapes by making some a^ixopriate estimate as to an average radius. 

B. Survey of Radiation x^gssure over a Broad Spectrum of Wavelengths 

In Table III-l, we have characterized the size of acceleration achievable from 
a broad spectrum of electromagnetic radiation sources. To do this, we have made 
use of the fact that maximum acceleration is achieved with Of = 1 

and for this value of OC , the efficiency factor for strongly refracting and 
absorbing materials of the order of unity. The power density 

was chosen as U = 1 watt'cu^ and the particle density as ^ ss / gram/cui^ 


Table III-l 


Frequency Range 
or 


Nominal 


Min. Part. 
Dia. for 
Good 


Acceleration in 
a vacuum of the 
min. particle in 
a 1 watt/cm^ 
radiation field 



' Microwave, X band 

8 - 12.5 GHz 

2 cm 

0.75 tfA) 

0.00067 cm/sec 


-Microwave, K band 

18 - 26.5 GHz 

1 cm 

0.25 cm 

0.002 cm/sec2 

9 

Microwave, 0 band 

60 - 120 GHz 

. 3 cm 

0.075 cm 

2 

0.0067 cm/sec 


* Infrared, submil- 
limeter 

1.2 x,l0^^ Hz to 
10l3 Hz 

2500,4/ 
to '3^ 

630^^ 

.0079 cm/sec2 
to .67 cm/sec2 

— 

Infrared, CO^ laser 

2.83 X 10^^ Hz 

10. 6 A. 

2.6^ 

1.9 cm/sec? 


Visible light 

6 X 10^^ Hz 

0.5/C 

O.l?^ 

AO cn/sec j 


Ultra Violet 

6 X 10^5 Hz 

o 

• 

o 

O.OlgA^ 

AOO cm/sec2 j 


* X-rays and rays 

10^7 to 10^2 Hz 

30 X to 

1.0 9 

2 * 

.0026 cm/sec { 

1 


0.0003 9, 

(atoms) 

1 


Note that the power density in these regions which is currently available from state- 
of-the-art power sources is generally —uch less than I watt/cm so that in practical 
equipment the accelerations achievable in the asterisked rows would be proportionally 
less than the values shown. X-ray pever is about 10"2 watts/cm as shown on page 8, 
and except for a few experimental lasers at a few hundred microns, there are no good 
sources in the submillimeter region. 
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It is evident from ti.e table that even the highest inlcrowave radar frequencies 
can in general produce only very small accelerations on the smallest particles 
with which they can interact efficiently. The acceleration values can be 
bettered somewhat if a higher intensity than the nominal C.W. power intensity of 
1 watt/cm^ can be achieved either by use of higher power sources or by focussing 
of the energy. 

The submillimeter infrared region might produce usable accelerations if 
reasonable power levels were available, but the means of generating and guiding 
submillimeter waves are very limited, and the power levels available decrease 
rapidly as wavelength decreases in this region. 

It would seem that because of their short wavelength X-rays would be a good 
source of accelerating individt 1 atoms and molecules because of the favorable 
value of the parameter x. However, the high energy of the X-rays makes them 
very penetrating and consequently they do not interact strongly with Individual 
atoms and molecules. The lead atom, which has relatively a very large X-ray 
cross section, has a cross section for emission of photoelectrons (the largest) 
of only 2,7 x 10“20 cm2 ^ for 10 KeV X-rays 1,24 8), This corresponds 

to an efficiency factor of only Qpr-= 10"^, Other atoms with smaller atomic 
numbers have correspondingly smallei effective cross sectiorfi. For example,* the 
air molecule has a cross section two orders of magnitude smaller than lead.^ 
Furthermore, the available power levels for X-ray sources are also very small. 

The 1 watt/cm^ value used to compute the accelerations in the above table is 
probably at least two orders of magnitude larger intensity than can be con- 
veniently achieved. For example, the Coolldge tube, long the standard source 
of X-rays, produces only 3.9 x ID'S watts/cm2 at a distance of one meter or 
3.9 X 10“3 watts/cm2 at 10 cm from the target with an accelerating potential 
of 100 keV and a current of 10 ma.& Therefore, the acceleration of 
0.026 cm/sec^ in the table represents a value that probably could never be 
achieved. In addition, inspection of Fig. II-l reveals that based on thermal 
energy arguments alone the minimum acceleration needed for gas molecules is approxi- 
mately 10^ cm/sec^. Consequently, we are well Justified in not pursuing the 
Issue of x-ray sources further. 


'C. Radiation Pressure from Sources in the Near-Ultraviolet. Visible & Infrared Rseior. 
Non Absorptive Materials 

Having focussed attention on the range of wavelengths near the visible, we will 
now show the results of calculations of particle acceleration due to radiation 
pressure as a function of particle size, refractive index, and wavelength of 
the source radiation. For each set of particle parameters considered, the 
calculations were done for three different wavelengths of source radiation, each 
corresponding to a known laser frequency; (1) .32^<(,(HeCd) , (2) .6328^41, (HeNe) , 

(3) lO.b^^^CCOp . The calculations were made in two parts, the first correspond- 
ing to nonabsorbing particles (real refractive index) and the second to absorbing 
particles (complex refractive Index) . 
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As noted earlier, the process of calculating the acceleration due to radiation 
pressure amounts to the determination of the efficiency factor for radiation 
pressure (j?pr ♦ In general, for a sphere, this amounts to solving Maxwell's 
equations where one expresses the incoming wave, the wave inside of thf 
and the outgoing wave In terms of an infinite sum of spherical Bessel ^ij-ictior 
and then matches boundary conditions at the surface of the sphere. TIjas was 
first done by Gustav Mie in 1908.7 For the work to be presented here, <;?P^^was 
obtained by piecing together and convolving separate more restricted theories 
of scattering such as the theory of Eeyleigh scattering® and anomalous 
diffraction^ with scries expansions and numerical calculations^^ derived 
directly from the formal Mie scattering theory. 


The initial acceleration as a function of rarticle size for non-absorptive 
spheres with refractive indices of m = 1.1, 1.3, 1*5, 1.8 and 2 for the 
three above noted wavelengths is shown in Figs. 111-1,-2, and -3. The particle 
density is ^ 1 gram/cm®, and the radiation Intensity is W = 1 watt/cm^. The 

dashed lines in these figures are a reproduction of the requirement which is 
shown in Fig. II-l for the minimum initial acceleration necessary to confine 
approximately 70 percent of the particles (i.e. 1-^ A*2.7/?)in a distance of 1 cm. 
It should also be noted that the acceleration is plotted in units of the earth's 
gravity ( G * 9.80 x 10^ cm/sec^) * 

The most Interesting feature of the results shown is the peak in the particle 
size dependence of the initial acceleration. For a given wavelength,* the peak 
is approximately independent of refractive index but diminishes in size as m-1 
approaches zero. Tlie position of the peak occurs for a particle radius that 
satisfies the requirement that /X « ®^^=Vand the maximum acceleration is 
approximately given by the following empirical formula 

where ^ in this expression is in microns and 

where the refractive index is constrained to be in the range t ^ |V\ jL 
The particle radii for which the acceleration is ^.707 of the maximum value 
(given above) is approximately given by the following empirical formula 


where 


( < a. in-5 


The accuracy of the above two formulae is /\/10 nercent. For particle radii 
such that X4.0.5 (eq. for * .3254{, a <.02^; for X “ •§32^, a 4, .0^; 
for 10.64<, a <^0.8^) the acceleration falls off as a^. The reasen for 
this is that for these small values of x the scattering Is essentially 
Rayleigh scattering® so that (pp^has a dependence of a^ but the 1/a in Eqn. IZX-3 
reduces this to the a^ dependence mentioned above. For particle radii such that 
X ;> 4 (eq. for .325^4, a 0.2^;for/\« .632^, a ^ .^)for 10.^. 
*^^1/0 acceleration falls c-£ as l/a because for these larger values of x 
Qpr is approximately constant' so that the 1/a dependence in Eqn. IZZ-3 prevails. 
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The range of reft i ivr indices f is sufficient to encompass almost 

all non-absoi'ptlvt "materials. The refractive index for most substances can 
be found in the Handbook of Cliemistry and Physics. Refractive indices for 
some representatli'j?. js:.'". * . i are given J.n Table III-2. 


Table III- 2 


Ice 

m= 1.31 


Uatcr 

me 1.33 


Sulphur 

me 1.99 


Chlorine (liquid) 

me 1.385 


Hydrogen (liquid) 

me 1. 11 


Aluminum oxide • J:., . 

me 1.67 


Sodium Iodide 

m- 1.77 


Sodium Chloride ' 

me 1.53 


Albite glass 

me 1.49 


Amber 

me 1.55 


Ivory 

m- 1.54 


Quartz 

m» 1.5 


Virus 

me 1.6 

(a - 005^ ) 

Protein 

m« 1.6 

(a e O.yc) 

Bacteria 

m« 1,4 

(a e l.C^ ) 


Although the refracti?j 2 irdex of most non-absorptive materials is not a drastic 
function of wavelengtK.in the region we are dealing with, there usually is some 
change ( /v«30%) wit.v.'..-3velength, particularly in the infrared. .The dependence 

of the refractive index on wavelength for a number of materials can be found 
in the American Institute of Physics Handbook. 
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Absorptive Materials 

The task of obtaining the efficiency factor for radiation pressure is 

more difficult for absorptive materials because^ except for small values of 
the parameter , there is no recourse other than the formal Mie theory. 

However, the task of carrying out a calculation using this complex theory is 
beyond the scope of this study. Consequently for the X=. 0.3;!^ _ 

and ^ <1? • 1 32^^ v;ave lengths we have only treated one refractive index^ 

1 .2,7 — 1*37 iL where i = • This refractive index is actually that for 

iron at "X” .H*j^ 4 ^and is one for which (ppr has been calculated.!^ It, 
however, is not precisely equal to the refractive index for iron at the specific 
wavelengths considered here ( eq. for X” « 3*2.5/^^ ; 

and for /77 .= /. 7- / 8iC\ 12.15 

but is typical of the refractive indices or many metals in this range of 

wavelengths, and thus serves as an illustration of the type of behavior to 
be expected. With this introduction, the initial acceleration as a function 
of particle size for the wavelength X- *3^1^.^ and shown in 

Fig. IIIr.4. The dotted line in the figure is the initial acceleration needed 
to achieve a 1 cm concentration density of the type described earlier, but 
for this case the temperature is taken at lOOO^K* This was done because 
lOOOOK appears to be a more realistic operating temperature than 300°^ since 
the materials are absorptive and will heat up when illuminated with high ■ 
intensity radiation. The initial acceleration as it appears in the figure 
has been normalized with respect to the density of the particle material so 
that the figure can be readily applied to a variety of materials. Tp find the 
initial acceleration for a specific material, one must only divide the acceleration 
appearing in the figure by the specific gravity cf the material. 


In common with the previous results for non-absorptive spheres, there is a 
peak in the acceleration dependence for values of particle size where the 
parameter oc “ — / • Similarly, at larger particle radii, the 

acceleration falls off as 1/a since the efficiency factor for absorptive 
particles also becomes constant. Th^ main difference between the absorptive 
and non-absorptive particles is that the acceleration does not fall off at 
small values of particle radius but approaches a constant value. This is due 
to. the fact that while the efficiency factor for scattering diminishes as a^ 

for small particles (i.e. f)/( ^ »S )> the efficiency factor for aoserption dimi- 
nishes only as af®^ 

The situation changes quite dramatically in the Infrared region since the refractive 
index of metals becomes quite large and is dominated by the conductivity of 
the material. It is given for ^ as^® 


y 

('30«r?v) (l-i) 
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where O’ is the conductivity in MKS units and . WhenYYVOC^ 10 , 

the scattering dominates the absorption, and the spheres can be treated as 
perfectly reflecting for purposes of calculating the radiation pressure. 

* The complex part of the refractive index is sometimes referred to as the 
absorption index, and reflects the fact that the material is conductive and 
thus can absorb energy. 
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For this case, all metals (where rA»jf> 10) will behave identically. This require- 
ment that n\' 3 C *>10 is essentially a requirement that the skin depth be small in 
comparison to the particle size. For the case of|MY'<,l absorption dominates the 
scattering and the theory becomes more complicated with different metals behaving 
differently depending on their conductivity. The initial acceleration for ^h^ 
case of iron (but normalized with respect to density) at shown in 

Fig. III-5. As can be seen, there are two regions of interest, one dominated by 
scattering and the other by absorption. In the scattering region, the normalized 
acceleration shown would be Identical to that for metals which have a conductivity 
at least as large as iron, since the efficiency factor for this portion of the 
curve was taken from that of perfectly reflecting spheres. The position of^ the 
peak in the scattering region also occurs for particle radii where A 

The region of the curve marked absorption is much more complicated and is unique 
to iron. In this region, the particle size is comparable or smaller than the 
skin depth. A complete discussion of this region is found in Van De Hulst. 

An additional constraint is Imposed in the case of absorptive particles, in that 
there is a limit to the amount of Incoming radiation intensity that the particles 
can accept in order to avoid melting. We have calculated this power for iron 
spheres on the basis of equating the power absorbed from the radiation to the power 
lost to black body radiation^^ at the melting temperature of iron That is 


Power lost 


- art 

Power absorbed =* •'W 


m-7 


III-8 


where XfcT^ V is the Stefan-Boltzmann constant, T is the terrperat 

the efficiency factor for emission at the center of the range ot 
wavelengths corresponding to black body radiation at temperature T, is 

the efficiency factor for absorption at the wavelength of the incident radiation, 
a is the particle radius^ and W is the intensity of the incident radiation# In 
calculating the radiant power lost to black body radiation, Kirchhoff’s law 20 
states that the efficiency factor for emission (the emissivity) at a particular 
wavelength is equal to the efficient factor for absorption at the same wavelength. 
This is the reason that the factor included in Eqn. III-7 for the 

emissivity. With this explanation, the maximum incident radiation intensity in 
order to avoid melting iron spheres is plotted as a function of particle radius 
in Fig. III-6 for .4S^4(and in Fig. III-7 for . 

These plots are superimposed on plots of the initial acceleration with the original 
power level of 1 watt/cm^ taken from Fig. III-4 and Fig. 111*5. 

The dependence of the maximum power on particle size is due to the dependencies 
of (?<i 4 (>i^and As can be seen, the initial acceleration can be greatly 

increased in all cases before melting will occur. The allowable power level is 
so high in the case of 10. 6 radiation that it is unlikely it will be approached 
In practice. The large size is due in part to the large size of 
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the refractive index at lO.^/ynaking the iron spheres seem highly reflecting^ 
but it is also due to the fact that the black body radiation occurs at a much 
shorter wavelength ( > which places it in a region where the efficiency 
factor for absorption, and nence emission, is much larger. The size of the 
allowable power levels (and the fact that high power levels can be achieved 
with the CO^ laser) and the lower required initial accelerations for large, 
particles (Sec Fig. II-l) more than compensate for the smaller accelerations 
obtained with the calculated power level of i watts/cm^ seen in Fig. IH-s' ' 
and ••ig. Ill- 3. 


Resonant Radiation Pressure 

The radj...::ion pressure on atoms and molecules was also considered. However, 
the light scattering formalism as treated above for pirticles in the micron 
size range results in accelerations much too small when compared to thermal 
energies from atoms and molecules. Referring to Fig. II-l, it can be seen that 
to satisfy the requirement based on thermal energy (curve 1) for a particle with 
a radius of 3 X (atomic size) would require an initial acceleration of«^105 g's. 

As will be seen below, accelerations on this order can be obtained from atomic 
resonant absorption in the visible wavelength range. 

21 

The force resulting from resonant absorption is given as 


F 



III- 9 


where h is Flancks constant, X is the wavelength corresponding to resonant 
absorption, ^ is the natural lifetime of the excited state, and f is the 
fraction of time the atom spends in the upper (excited) state. The quantity f 
is linearly dependent on radiation intensity at very low power levels but 
eventually stops Increasing with increasing power since ultimately the atom 
can only absorb a photon until after returning to the ground state, a time delay 
which is equal to the lifetime of the excited state. At the point when the 
absorption does not increase with power level, the line is said to be saturated. 

This equation has a simple interpretation. The quantity •n/X is the mumentum 
of a photon and is the number of photons per second tnat the atom can 

absorb as limited by the lifetime of the excited state and the 

amount of time the atom spends in the ground state without absorbing a photon 
(1 - f ) . The size of this fon-.e can be estimated by letting 

(the value for the sodium D line), and f » 1/2 which assumes the resonance line 
Is saturated and the ground and excited states are equally degenerate. 

The resultant force is 6.6 x 10"^^ dynes and the initial acceleration for the 
case of the sodium atom v;ould be 1.76 x 10^ G's. The power level to produce 
saturation in this case is .019 W/cnr (as will be shown below) . 

Equation III- 9 can be written in an alternate form by noting that the natural 
linewldth of the resonance line is related to the lifetime of the excited state, 
i.e. ^ l/sittll . If we define a quality factor for the resonance line 
as ^ sVo/A)) > rewrite the force equation as 
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vhere c is the specr'j pf If we assume that the quality factor of most 

resonance lines is aj'AV ' ’ independent of the resonance frequency, ^ 

this equation makes it clear that the force resulting from atomic resonant 
absorption is proportional to the second power of the resonant frequency- Con- 
sequently, one should not expect significant forces (or accelerations) from 
resonant absorption in the microwave frequency range* 

21 

Ashkin makes the observation Jn his paper that, depending on the applica- 
tion, doppler shifts of the atoms as they gain speed may limit the achievable 
acceleration by causing t*he atoms to move off of resonance* One solution 
to this, as Askin xs to power broaden the resonance line so that the 

effective linewidth i»^%inc recced. A broader resonance line may also make 
frequency control of t:he laser supplying the resonant radiation unnecessary* 

It then becomes relevant ‘to know how much power is needed to achieve saturation 
since broadening ensues at powers greater than this. To do this calculation 
we need to know the scattering cross section for resonance absorption. The 
optical theorem states 


Cu= 
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vhere is thi:'.«ttu:il cross section for absorption, k is the wavevector 

of the incoming radiation 4 Pi stands for^he real part ofV^s/Tj' , 

and is the polarizability. Using a simple classical model for damped 

harmonic motion of an electron^the polarizability is given as 23 
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vhere is the resonance frequency, M is the electron mass, JL is the electron 

charge,and is tb.i U«.„rfidth. The cross section for absorption is thus 


- - „i.i3 

Evaluating this on resor.o:ice by taking yields Ce.ts = 

ItOGX The force on the atom based on the absorption cross section is given 


/O 
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Equating this equation to Eqn. III-9 and then solving for the intensity W 
gives W-\/ /,9X/£rc/yi/>i^ Of' ‘Ot°t /»!*■ to saturate 

the line. This intensity i? sufficiently low that saturation can easily be 
achieved and significant power broadening can be accomplished with realizable 
laser sources. Ashkin discusses the size of the broadening as a £unction^(^ 
power. level in his paper. 



The real limitation to the use of resonant radiation pressure lies in the 
requirement that a laser source be available at precisely the resonant 
frequency of interest. However, strides are being made in the development 
of tunable lasers so thac this interesting technique may lend itself to 
applications in the near future. 


Conclusion 


This completes the section on radiation pressure. However, since the size 
of the acceleration that can be produced is directly related to the intesity 
of the "light" source some comment as to the power levels that can be achieved 
is in order. For the ultraviolet Helium Cadmium laser (.32^) only 10 mW 
total power is possible. The Helium Neon laser is capable ^ 100 mW^^ Xhe 
C0_ laser is capable of extremely high powers, 10 kW or greater, but these are 
only for very expensive, very large lasers. Considering a laser 3 ftet long 
as a maximum size, 100 W Is a reasonable power level that can be expected. 
All of these values are for CW operation. It should be noted that very high 
intensities can be achieved by focussing the laser and thus placing the total 
power available in a small area. The theoretical limit to this focusing Is 
an area that Is equal to the square of the wavelength. 







IV. SOUND PRESSURE 


Our investigation has been United to conpressional waves in gases and 
liquids, and the effects of such waves on particles suspended in the nedlurn. 

Particles nay be moved by two effects: radiation pressure and/or acoust\c 

streaming. 


Radiation Pressure 


The DC radiation pressure is described by an equation similar in form to the 
expression for radiation pressure in an electromagnetic wave, that is, the 
pressure is equal in value to the energy density of the wave.^^ The net DC 
force produced on an obstacle in the medium carrying the sound wave is equal 
to the effective cross section area ratio for momentum transfer times the 
actual cross section area of the particle times the energy density of the 
sonic wave. As in electromagnetic radiation the effective cross section 
area of the particle is much smaller than the actual cross section for par- 
ticle sizes smaller than a wavelength, so that little force results in that 
regime. 

Now in gases, the maximum sonic frequency which can be transmitted Without 
excessive attenuation due to the medium itself is in the order of 10^ Hz, 
corresponding to a wavelength of 10"2 to 10"1 cm.27 Thus, manipulation of 
small particles in a gaseous medium using sonic radiation pressure appears 
to be impractical, since micron sized particles would be much smaller than 
a wavelength of sound. 

9 

In a liquid medium, however, frequencies up to the order of 10 Hz can ue 
propagated, corresponding to wavelengths of 1 to 10 microns, so that smaller 
particles can be moved in a liquid medium. Ne have not been able to pursue 
this use of sound radiation pressure in dutail, because of limits of time and 
money, but it is worthy of 'further investigation, especially as to the char- 
acteristics of practical sources of energy at such frequencies, and the sonic 
reflective or absorptive characteristics of particular substances. 


Acoustic Streaming 

Acoustic streaming is another sonic effect which we have not been able to 
pursue in detail during this study, but which shows Interesting possibilities 
for manipulation of small particles. Acoustic streaming is a DC, convective 
type of motion of the medium which occurs under certain circumstances with 
sonic radiation. Particles suspended in the medium will experience degrees 
of motion related to their size, thus forming a possible means of sorting. 

One type of acoustic streaming is caused by the presence of solid boundaries 
such as the walls in a resonant chamber containing a standing wave. 

The streaming or vortical motion of the medium results from the fact that no 
slip motion can occur at the walls. A well known example of this phenomenon 



Is the *'4^ . *• 

shaken into a 
tube is. -oner ^ 
standing 

and deposited. out 
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'fine sand or lycopodium powder is 
• stVpoing wave visible. When the 
Icivrj from the interaction of the 
•c, i,yJ)e,.V^sjiises the powder to be transported 
of th^ standing wave. 


Lord Rayleigh and among otViers, solved this problem for tKe 

particular case of a standing; wive of. sound existing in a waveguide formed 
by a pair of parallel plffifc?'1t?rrili;ated at one end by a reflecting wall so 
as to produce a standing wave.^ They derived equations for the net horizon- 
tal and vertical DC moticns of the particles of the medium during 

steady state excitatie.-.. 




Outside the motionless thin film immediately adjacent to the walls, the net 
horizontal velocity, r.-- -i" given by 


u„ 


3u ^ sin 2 Kx 
o 

16c 



( 1 ) 


and the net DC vertical velocity, perpendicular to the walls, is given by 


-3u 2K v s- 2Kx f 
o I 


(71 - y>” 


16 a 




•} 


( 2 ) 


where 

y Is the distance from one wall to the point under consideration 

y. is the distance from one wall to the mid-plane of symmetry 

Is the maximum value of the instantaneous horizontal velocity 
K° is the wave number of the sound wave, i.e, K « (2jt/X) 

X is the horizontal distance from the sound source 

c is the velocity of sound in the medium 


It is seen that these motions describe a series of vortices periodic with 
respect to x in the distance X/2, The vertical velocity, V 2 , is negative, 
that is, towards the walls, at the loops, and positive, or away from the 
wall, at the nodes, vanishing in each case both at the wall, y « 0, and at 
the plane of symmetry, y « yj^. The equation for U 2 > the average horizontal 
velocity, shows that it is zero at both the nodes and the loops, and has 
maximum positive and negative values half-way between the nodes and the 
loops. 

The nature of the vortices explains the behavior of the Kundt's tube. It 
is seen that particles suspended in the medium will be carried horizontally 
to, and deposited out at, the nodps. 
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In Figure 1, AB represents the wall and CD the plane of syranetry, and the 
directions of motions in the vortices are indicated by the arrows. 



Figure 1 - Streaming Vortices in a Chamber Containing a Standing Wave 


We have in a preliminary way extended Rayleigh* s analysis to explore the 
behavior of the vortices in the frequency domain and to determine the ef- 
fect of the streaming motions of the medium on particles of various sizes 
Immersed In the medium. 

One can assume a transient situation in which the sound is turned on and a 
short time later the streaming motions of the medium reach steady state 
velocitit . At this early time the suspended particles, which are assumed 
to be very much larger than the particles of the medium, are assumed to 
have reached only a small fraction of the velocity of the medium. We can 
then calculate the force on these particles using Stoke* s law which states 
the proportionality between drag force and the product of the differential 
veloc*ty, viscosity of the medium, and particle radius. Dividing the drag 
force by the mass of the particle yields a value for initial acceleration 
of the particles. 

Initial accelerations computed in this way do not of course show directly 
what can be achieved in a particle sorting device based on these principles, 
but the order of magnitude of the initial acceleration produced and its var- 
iation with particle parameters gives a general indication of the potential- 
ities of such devices. 

The results of this preliminary analyses yield the following conclusions: 
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1) At low frequencies of excitation the particles would tend to be deposited 
at the nodes in relatively narrow bands« 

2) At high frequencies the deposition would still concentrate at the nodes 
but much more diffusely than at low frequencies. 

^ % 

r 

3 } The initial accelerations are greater at the higher frequencies* 

4) The acceleration produced is Inversely proportional to the product of 
the density of the particle and the square of its radius. 

5) In a practical case we assumed water droplets of one micron radius in air 
in a parallel plate channel containing a standing wave with acoustic sound 
level of 120 dB above 0.002 dynes/cm^. At low frequencies the maximum 
Initial axial acceleration was calculated to be 17 cm/sec^ or ~^.017 G, 
and the maximum initial vertical acceleration was found to he 39 cm/sec^, 
or 0.039 G. At high sonic frequencies the maximum initial axial acceler- 
ation remains about the same^ Imt the vertical components increase with 
frequency. 

The preliminary results indicate that acoustic streaming will produce signi- 
ficant motion of particles. Further work is needed to confirm what has 
already been done, and to work out in more detail the particle orations due 
to streaming forces in a viscous medium. 


Other Acoustic Forces 

30 

Three other acoustic forces are listed by Hueter and Bolt , but have not 
been investigated in this program. They are mentioned here as a guide to 
possible future work. They are as follows: 

Direction 


Type of Force 

Cause 

Drag Coefficient 

of Force 

"Average Stokes" 

Temperature 
dependence of 
viscosity 

3(y - 3)T, 

Toward 

Source 

"Oseen" 

Wave 

distortion 


Depends on 
phase 

"Bernoulli" 

Hydrodynamic 

Flow 

2 2 
.niHz 

Attraction 
between 2 
particles 


The symbols used are defined as follows: 


r ■> particle radius 

d B distance between two adjacent particles 
Y “ specific heat ratio 
Tj ■ viscosity coeficient 
p » density of the medium 
C 2 ■ 2nd harmonic content 

^ • phase shift between fundamental and 2nd harmonic of distorted sound wave 
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V. PARTICULATE RADIATION 


We consider particulate radiation to denote strccuns of elementary particles 
such as electrons# ions# or uncharged gas molecules. 

We have considered the case of an electron beam impinging on an aggregation of 
particles and have calculated the forces and accelerations produced by relatively 
low energy electrons whose mass is approximately equal to the rest mass.^ViCe 
will assume fcr simplicity that the electron is totally absorbed by the particle 
it strikes so that it gives up all its energy and momentum to the particle. 

The derivation of the force equation is quite simple, and it proceeds as follows: 


The force on the particle is just the momentum of each electron times the nusher 
of electrons captured per second 



Ij-Zi'xlO 

fer el&c'hroA ^ec Crtt** 



Airea <j-f- 


T-l 


For velocities 
given by 


much less than the velocity of light the electron velocity is 

« 

5 *. 92 . ycio'^ c^^sec, 


where V is the potential difference in volts through which the electron was 
accelerated. 
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The middle term arises from the fact that there are 6.28 x 10 electrons/ 
second per ampere of current. J is the current density in amocres/cm^. 

Putting in values for K^ao^i^we get for the force on a spherical particle 
Of radius a cm 


F = )M xla'^lfV' T-3 


The acceleration produced in a vacuum is obtained by dividing both sides of the 
eqviation by the particle mass, (4/3)Trd^^yielding 
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In a specific case we assur.e an electron bean of 100 volts energy and 10 
millianperes/cn2 current density, in 5 >inging on particles of 1/t radius and 
1 graiVcn^ density in an evacuated chairher. Putting these values into the 
acceleration equation, we find a large acceleration of 254 cn/sec^ or 0.259 G 
for this example. 


li 


Interaction Chamber 
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Ihe above calculations assumed an idealized situation in which electrons 
continuously strike particles and give up all their energy and monientun to them. 
In a practical case, however, there are several characteristics of the electron- 
particle collisions which must be considered before it can be determined whether 
the idealized model can be practically realized. Let us consider the case of 
the evacuated chamber first and determine what happens when the electron beam 
is turned on. We will assume em electron energy in the order of 10 volts, low 
enough that secondary electron emission will not complicate the problem. 

Initially the electrons will be absorbed by the particles pretty much in accord- 
ance with the simple model on which the force equation was derived. As time 
goes on the particles will become negatively chcurgad, and no further electrons 
will be absorbed when the potential of the particles reaches the energy of the 
beam. 

Nevertheless, even though no further electrons are absorbed by the particles 
after charge equilibrium is reached, additional electrons will be shopped or 
scattered by the potential field of the particles, so that momentum should 
continue to be transferred to them, and the effective cross section of the 
peurticles should continue to be the same order of magnitude as their actual 
cross section area. But then the problem becomes more complicated because 
potential differences may exist between the walls of the container and the 
particles which will produce additional forces on the particles. 

It was beyond the scope of this study to work out the kinds of interaction 
chambers which might be required to utilize the force created by an electron 
beam inpinging on particles. The charges produced on the particles certainly do 
complicate the practical problem of utilization devices compared to the use of 
the other force fields considered, which do not produce charges on the p^u: tides. 

It should also be noted that it may be difficult to use an electron beam to 
numipulate particles in a viscous medium such as air because of the large 
absorption or scattering of the beam by the gas molecules themselves. For 
example, tlie attenuation constant for a 100 volt electron beam entering neon 
gas (which has a relatively low collison cross section for electrons) at 1 torr 
pressure is 10 per cm, so that the beam would be reduced to 1/e of its initial 
intensity in 0.1 cm, 31 
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VI. ELECTROSTATIC / ‘ 


Electrostatic fore. . *t ‘ t of coulor.b forces , i.e., forces on charged 
particles in an elc\:r:^^ or electrophoretic forces , i.e., forces on 

small charged parti J * > Vodial solutions, or dielcctrophoretlc forces , 

i.e., forces on uneb/.* *. • teles in an electric field. ^ V 

We have not treated coulomb forces, because in general, they are large forces, 
their characteristics are widely known and utilized in electrostatic pre- 
cipitators, and they are not generally capable of fine control because of the 
uncertainty of the amoCiU charge put on the particles by corona discharge. 
We have not considered electxophoretlc forces because we understand that NASA 
already has a program of evaluation of electrophoretic effects in the zero 
gravity environment. ^ - 

Dielectrophoresis *** 


Dielectrophoresis is theT^rae given to the production of forces on uncharged 
particles in a non-uniform electric field resulting from the electric po*lari- 
zation of the particle. The polarization is caused by shifts in the electron 
clouds of the atoms produced by the electric field so that each atom acquires a 
dipole moment. If the field is uniform, equal and opposite forces ^are produced 
on the induced charges of the dipole, and no net force results on Che particle. 
If the field is non-uniform, however, the forces on the dipole are not equal 
and opposite, and a weak r. t force exists on the particle. Dielectrophoretic 
forces are predictable an.; . reproducible, since they originate from fundamental 
properties of the particle and of the field. 


The derivation of the expression for the magnitude of the dielectrophoretic 
force is given in advanced texts on electromagnetic theory, such as Stratton's 
Electromagnetic Theory, 1st edition 1941 published by McGraw Hill. The deriva- 
tion will be sketched out here, and then some typical examples of the magnitude 
of acceleration which can be produced on various substances in a simple coaxial 
structure with a potential difference applied between the coaxial elements will 
be given. 


Stratton gives the change in energy resulting from introducing a body of per- 
mittivity ^2 a uniform medium of permittivity as 

ir~i 




where E is Che electric field vector inside the body, and Eg is the electric 
field vector in the neighborhood of the particle before it was brought into 
the field, i.e., Che applied field. If the body is spherical, Stratton shows 
the relation between ^ and Eg for a uniform field to be33 
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Looking ahead, one notes that the field must be non-uniforin in order for 
there to be a force, so that there is a question as to the validity of 
using Eq. VI-2 for this case. It can be shown, however, that Eq. VI*2 
gives the zero order tern in the electric field for the non-uniform case^ 
and also becomes a reasonable approximation to the total field if the 
variation in the field is small over the volume of the sphere. If we 
now substitute the value of E given by Eq. VI-2 into the integral ex- 
pression for the energy, given by Eq. VI- 1, v;e note that the integrand 
is spatially varying, thus complicating the evaluation of the integral. 
Nevertheless, the integral can be evaluated by taking the electric field 
to be constant, using the value at the center of the sphere as the mean 
value. This will be a good appreximation to the integral if the varia- 
tion in field is small over the volume of the sphere. i 

The resulting value of the change in energy due to introducing the sphere 
turns out to be 


U 






IP- 3 


where a is the radius of the sphere. 


The force acting on the sphere is the negative of the gradient of the energy 
change in the direction of the force. T! us 


F--VU 


IPP-^ 


and performing the indicated differentiation, and substituting for£| 

and ^ 2^0 ^l ^2 relative dielectric constants of 

the medium and the particle respectively, and €, is the permittivity of a 
vacuum, we find that the expression for the force on the particle is 


p = S.-vS. ET - £ 


This expression is in HKS units. The force is in newtons (1 newton = 10^ 
dynes); the particle radius a is in meters, the electric field, E^, is in 
volts per meter, and is the permittivity of a vacuum = 8.855 x--10"^^ 
farads /me ter. 

The main assumptions on which the validity of the above analysis dr-pends are 
that (1), the particles are small compared to the distance in which the elec- 
tric field changes appreciably, and that (2), the number of particles per 
unit volume is small enough that the presence of the particles does not 
appreciably affect the pattern of the applied field. For small particles 
dispersed in a reasonably sized interaction chamber, these assumptions 
should be met . 







Exatnples of Dielectrophoretic Force 

Let us consider the case of a particle**fleld interaction chamber consioting 
of two coaxial electrodes of radius and R 2 meters with the particles 
located in the annular space between the electrodes, and a potential dif-^ 
ference volts applied to the electrodes. Such a configuration yie^t^s 

a highly non-uniform radial electric field between the electrodes. It can' 
be shown that in this configuration the product of where )r\ is in 

the radial direction in this case^ is 

£ v£, 

V Ctiii - , _ - , —,2. 

Zf/Z/f J 

where n is the radial distance outward from the axis of the cylinders. 
Thus in this particular configuration the force on the particle is given by: 


-^TTa 


1^-7 




and the acceleration in a vacuum by * 

/i , I — • 

Hccel^- ^ 
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where ^ is the density of the particle in Kg/meter'*. 


Note that the acceleration in a vacuum produced by the dielectrophoretic force 
is dependent on the dielectric constant of Che particle and its density, but 
is independent of its size. 

If we now assume that the applied potential, is 10,000 volts, that R]i is 
1 cm and R 2 is 2 cm, that x = Rj^ = 1 cm (meaning that the particle is located 
near the surface of the inner electrode), we find accelerations expressed in 
units of g as follows for a number of substances: 



I 


1 




lable VI- 1 
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Oielectric 
. f: * . lant^ 

(dimension- 

less) 

1 (K,-!) 

Acceleration produced ai 
the center of a coa>(ial 
chamber with Rj = 1 cm, 
R 2 = 2 cm, and 10,000 
volts applied 

Substance 

• . t 

i 

Density, 

in , 

Kgr/meter-' 

0^2+2) 

Clausius** 

Mosotti 

Constant 

Water 

1.0 X 10^ - 

81.0 

0.97 X 10'^ 

-0.55 g 

Oil 

0.8 X 10^ 

2.2 

0.36 X 10"^ 

-0.21 g 

Paper 

~ 0.8 X 10^ . 

^ J 

■ -2.5 

0.42 X 10"^ 

-0.24 g 

Boro- 

Sllicate 

Glass 

2.24 X 10^ 

4.86 

0.25 X 10”^ 

-0.144 g 

Diamond 

3.51 X 10^ 

16.5 

0.24 X 10’^ 

-0.135 8 

Mica 

3.0 X 10^ 

6.0 

0.21 X 10"^ 

-0.119 g 

Sulphur 

2.07 X 10^ 

4.2 

0.25 X 10"^ 

-0.142 g 

Dry Wood 

■V 0.5 X 10^ 

~ 6.0 

1.25 X 10"^ 

-0.712 g* 


If we assume that a lower limit of manipulation would occur when the energy 
added to the particle by the dlelectrophoretlc force acting on the particle 
all the way from the outer electrode to the Inner electrode is equal to its 
thermal energy, we find an estimate of the minimum size of manipulatable 
particles to be about 0.06 microns for the particular interaction chamber 
described above. 

Suimnary of Dielectrophoresis Effects 

In summary, the accelerations produced by dielectrophoresis are very substan- 
tial, being in the order of 0.1 to 1 g, they are independent of particle size, 
they are inversely proportional to the density of the particle, and they are 
a function of the dielectric constant of the particle. Particles down to about 
0.05 microns radius should be manipulatable before random thermal energy inter- 
feres. 
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VII. MAGNETIC FORCES 


We arc all familiar with the very large forces exerted on ferromagnetir sub- 
stances by magnetic fields generated either by electromagnets or by permanent 
magnets. These forces nay be exerted on particles of ferromagnetic substances 
which have no initial magnetic poles as well as on permanent magnets. In the 
former case, the magnetic field causes the highly polarizable domains oi jthe 
magnetic material to line up with the field. If the magnetic field is uniform 
In the direction of the field, equal and opposite forces are applied to the 
induced magnetic poles in the particle, and no net force on the particle 
results. 

If, on the other hand, the magnetic field is non-uniform in the direction of 
the field, then the forces on the induced poles are not equal, and there is 
a net force on the particle which urges it in the direction of the stronger 
field. The direction and magnitude of the force are independent of the po- 
larity of the field, so that a net force is produced in AC as well as in DC 
fields . 

The derivation of the force equation for materials in a magnetic field fol- 
lows the same steps as the derivation for the dielectrophoretic force de- 
scribed in the previous section, and equations VI-1, 2, 3, 4, 5 and ^ when 
changed to magnetic units, describe the steps in the derivation of the 
magnt'tic force. , 

For the magnetic case, the dielectric constant, of the particle is re- 
placed by the magnetic permeability,^ of the particle, and the dielectric 
constant, of the medium, is replaced by unity, since it Is^assumed that 
the medium will be non-magnetic. The electric fields, § and E , are replaced 
by magnetic fields H and H . The resultant equation for the magnetic force 
on- a spherical particle is^ 

P =r vii-i 


where 

a 


is the magnetic permeability of the particle 
is the radius of the particle in cm 
is the field strength in oersteds 

is the directional derivative of the field in the direction of the 
field in oersteds /cm 


We note that this force is proportional to the volume of the particle, 
whereas electromagnetic radiation pressure, sound pressure, and electron 
beams all produced forces proportional to the area of the particles. 


Note that this is in cgs units, while dielectrophoretic force is in mks 
units. In the following, we will take the case where the field and its 
derivative are both in the x direction. 
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If we divtdg both sides of the equation by the mass of the particle, where 
M “ (4/3)j:a , we get the expression for the acceleration produced in a 
vacuum, or the initial acceleration produced in a viscous fluid. 


Accel 

X 


F 

X 

M 


0.238 


ICu » 1) 
p(p + 2) 


dH 

H ■ cm/ sec 

ox ax 


VII-2 




and we see that this is dependent only on the density and the magnetic prop- 
erties of the particle, and not on its size. 


Forces on Ferromagnetic Particles 

Kow if |i is much greater than 1, as is the case for ferromagnetic substances, 
the fraction Oi-l)/ (p+2) is practical unity, so that the acceleration equa- 
tion simplifies to . 


H dH 

Accel^ ** 0.238 — cm/ sec for ferromagnetic particles 


^VII-3 


It is possible in magnetic equipment to develop a product of (dH^^j/dx) 
of magnitude 107 oersteds^/cm, and putting this number along with a typical 
ferromagnetic material density of 8 grams/cm^ into the acceleration equation, 
we get 

2 

Accel (ferromagnetic) = 300,000 cm/sec VII-4 

« 300 G 

This is a very large effect indeed, and can be used, for example, to remove 
ferromagnetic particles from other particles with a good separation effi- 
ciency. 

Its magnitude is sufficient to warrant a look at the limiting velocities in 
viscous media such as air or water. To determine the limiting velocity, we 
equate the driving force as derived above to the drag force as given by 
Stoke* s Law. Again assuming a practical value of 10' oersteds^/cm for 
H (dH /dx) as in the previous example, we find 




6jt1|av 


driving force drag force 


v “ 5,3 X 10^ ~ cm/sec 
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where 7) is the vii' ii.y of the fluid in poises. For a ferromagnetic par- 
ticle of 1 micron radius the limiting velocities turn out to be 

(a-j's’ dt 30 cra/scc VII-7 » 

(watt/ - ■ ' cm/sec VII-8 ' 


These velocities are very substantial, even in water. In air they should 
produce usable .*;ffects even for magnetic particles smaller than 1 micron 
in radius. 





Forces on Diamagnetic and Paramagnetic Particles 

Most non- ferroriiagnetic materials have permeabilities which differ slightly 
from unity, typically by £» , ; *• parts per million. The permeability may be 
either slightly greater than unity, as in paramagnetic materials, or slightly 
less than unity, as in dl.anagnetic materials. 

« 

The force equation for such substances can be sinplified because the factor 
0i+2) in rhe denominator is very nearly equal to 3 for diamagnetic and para- 
magnetic substances. In this case the force equation becomes • 


H „ *^^ox 

3 ^ox dx 


dynes 


VII- 9 


and the acceleration equation for the practical case of H (dH /dx) « 10^ 
oersteds^/cm becomes 


Accel « r— «= 8 x 10^^ ■- cm/sec^ 

X M p 


VII- 10 


The factor (|j,-l)/p is called the "specific magnetic susceptibility," and 
values for this parameter for many substances are listed in the Handbook 
of Chemistry and Physics, 42nd Edition, pp, 1647 to 2657. Specific suscep- 
tibilities run from about + 10-4 to -10"o for a wide range of elements and 
compounds. Taking a few examples, we find values as follows: 


Substance 

Specific Magnetic 
Susceptibility 

Acceleration 1 reduced in a Practical 
Non-Linear Field where H (dH /dx)« 
10^ oersteds 

Mercury (liquid) 

-0,15 X 10”^ 

2 

-0.12 cm/sec 

Arsenic 

-0.31 X 10’^ 

-0,25 cm/sec^ 

Wood 

-0.3 to -0.7x10"^ 

-0,24 to -0,56 cm/sec^ 

Aluminum 

+0,6 X 10"^ 

+0,48 cm/sec^ 

Air at NTP 

+24 X 10‘® 

2 

+19 cm/sec (but this is too low a 
velocity to compete with random ther- 
mal velocity of the molecules) 
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In summary we can say that the forces produced on ferromagnetic particles 
by application of a non-uniform AC or DC magnetic field are very large. 
They will produce large accelerations In a vacuum which are Independent of 
the size of the particles. E«^en In viscous media such as air or water, 
appreciable values of limiting velocity can be obtained. 


'V 


The analysis gi\en here applies only to low frequency AC fields where the 
curl terms in Maxwell's equations are negligably small. This will generally 
be true for excitation frequencies in the audio frequency range. 

Although the forces on non- ferromagnetic materials are orders of magnitude 
less than the forces on ferromagnetic materials, they may be sufficient to 
produce useful sorting effects in an evacuated chamber. 
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VIII. SL1IMARY 


It was the purpose of this study to obtain quantitative understanding of a 
series of weak forces acting on small particles. The forces considered were 
those derived from light, heat, and microwaves (electromagnetic radiation 
pressure), electric interactions, magnetic interactions, particulate inter- 
actions, and sound. Throughout the body of this report, we have attempted 
to give a physical understanding of the forces considered as well as formulae 
which express how the size of the particular force depends on the physical 
and electrical properties of the particle. 

In the second section we have evaluated the drift velocity in a viscous fluid 
as a function of Initial acceleration in addition to considering the competing 
effects of thermal random motion. Out of all of this comes a means of selectively 
sorting or moving particles by choosing a force system and/or environment such 
that the particle of Interest reacts uniquely. 

The table below will serve as a summary of the forces considered and a demon- 
stration of how the Initial acceleration, drift velocity, and ultimate particle 
density distribution (/^ryCis the distance in which approximately 70 percent 
of the particles will be confined) is affected by particle, input, and environ- 
mental parameters. 

The dielectrophoretic and magnetic forces are identical as to their dependence 
on particle size and are unique in that the initial acceleration does not depend 
on particle size. The dependence of response parameters (acceleration, drift 
velocity, AOC) for the electromagnetic radiation force is probably the most 
unique of all the forces considered because of the strong dependence of the 
efficiency factor on particle size. Consequently, radiation pressure offers 
the best possibilities for selectively moving particles on the basis of particle 
size alone. Sound pressure is not included in the summary because the force is 
quite complex and not amenable to a short summary. In all of the cases in the 
table, the drift velocity and the parameter A'X are independent of particle 
density since the particle forces are independent of mass. 

It is difficult to summarize the sizes of the acceleration, drift velocity, and 
confinement distance because they are very dependent on the excitation levels 
used as well as the particle parameters. However, in order to convey some 
indication of the size of effect to be expected, values which are typical in a 
gross sense are given for the response parameters for each of the forces in the 
table. However, for complete appraisal of the sizes of the response parameters 
that can be achieved, particularly for the case cf radiation pressure, reference 
should be made to the appropriate section of this report. 

Tl^e reader should note that in the table the expressions for acceleration, 
drift velocity, and confinement distance under rev 2, Magnetic Force, are given 
In mks units. The equations used to describe n3g-..ctlc effects in the body of 
the text are given In cgs units. 
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^ = radiui of particle 


dezsizy of parciclc 
visccsity of fluid 


Boltzr i-r>ns constant 


= absolute teo<perature 
» permittivity of particle 
» permittivity of background 
= permeability of particle 
/. s permeability of free space 
So = electric field 
|4e** magnetic field 


s electric field directional derivative 


r = cagmetic field directional derivative 
efficiemcy factor for radiation pressure 




vavelemrth of E.M. radiation 


V/- intensity of E.M. radiation 


« speed of light 
= refractive index of particle 
“ Planck's constant 
= resonant frequency of atom 
= particle mass 

B quality factor of resonant line = 
B current density 
B accelerating voltage of electrons 
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